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S ■ Abstract 
OO , 

O ■ As a first step towards a nonperturbative investigation of tlie gauge-fixing (Rome) 

■ approacli to lattice cliiral gauge tlieories we study a U(l) model witli an action tliat 
^ ! includes a local gauge-fixing term and a mass counterterm for the gauge fields. The 

model is studied on the trivial orbit so that only the dynamics of the longitudinal 
gauge degrees of freedom is taken into account. The phase diagram of this higher- 
p • derivative scalar field theory is determined, both in the mean-field approximation 

D ! and numerically. The continuum limit of the model corresponds to a continuous 

■ phase transition between a ferromagnetic (FM) phase where the global U(l) sym- 
metry is broken, and a so-called helicoidal ferromagnetic (FMD) phase with broken 
U(l) symmetry and a nonvanishing condensate of the vector field. The global U(l) 
symmetry is restored in this continuum limit. We show that our data for the mag- 

53 ■ netization in the FM and FMD phases are in good agreement with perturbation 

theory. 



X 



Institute for Computational Physics, Humboldt University Berlin 
Invalidenstr. 110, 10099 Berlin, Germany 
e-mail: bock@linde.physik. hu- berlin. de 



Department of Physics, Washington University 
St. Louis, MO 63130, USA 
e-mail: 'maarten@aapje.wustl.edu 



School of Physics and Astronomy, 

Beverly and Raymond Sackler Faculty of Exact Sciences Tel- Aviv University 

Ramat Aviv 69978, Israel 

e-mail: ftshamir@wicc.weizmann. ac. il 



W. Bock, M. Golterman and Y. Shamir 



1 



1 Introduction 

All existing lattice fermion formulations have in common that they are in conflict 
with chiral gauge invariance. It is well known that for example the Wilson-term 
[|^], which is used to remove the 15 unwanted species doublers at the corners of the 
four dimensional Brillouin zone, is not invariant under chiral gauge transformations, 
because it has the structure of a mass term. 

Most fermion formulations can be rendered gauge invariant by inserting Higgs 
fields. The Wilson mass term for instance turns into a Wilson- Yukawa term which 
is invariant under chiral gauge transformations |]2|, |^ . These Higgs fields do not need 
to be added by hand. They appear automatically in the gauge noninvariant model 
when performing the integration over all gauge fields in the lattice path integral 
with the Haar measure Q. The Higgs fields can be identified with the longitudinal 
gauge degrees of freedom. 

The failure of most proposals for lattice chiral gauge theories, such as the Eichten- 
Preskill model the domain wall fermion model with waveguide ||^ or the Wilson- 
Yukawa (Smit-Swift) model is connected to the dynamics of these Higgs fields. 

In the Wilson- Yukawa model, for example, it has been shown, that the 15 un- 
wanted species doublers can indeed be removed from the spectrum by means of the 
Wilson- Yukawa term in the strong Wilson- Yukawa coupling symmetric phase. The 
model fails however because the left-handed fermion, which is the one that is sup- 
posed to couple to the gauge field, forms a fermionic bound state with the Higgs field 
that does not transform under the gauge group. This neutral left-handed fermion 
pairs up with the right-handed fermion to form a Dirac fermion which in the con- 
tinuum limit decouples from the gauge field 0. Later, arguments have been given 
that the spectrum in a symmetric phase is in general vector-like P, ^. 

The details of the mechanism which spoils the chiral nature of the fermions differ 
from model to model, but remarkably the Higgs fields (longitudinal gauge degrees of 
freedom) play the key role (for recent reviews see refs. ^, 0). It is therefore natural 
that one should try to use gauge fixing to control the effect of these longitudinal 



gauge degrees of freedom [[Tl]] 



Gauge fixing has been put forward some time ago as a method to discretize chiral 



gauge theories on the lattice |T^. It was proposed in ref. to use perturbation 
theory in the continuum as guideline and transcribe the gauge-fixed continuum path 
integral to the lattice. Since the fermion part in the action breaks chiral gauge 
invariance, the lattice model is not invariant under BRST symmetry. The symmetry 
breaking terms are compensated by counterterms that have to be added to the 
action. In four dimensions only the counterterms with dimension smaller than or 
equal to four need to be considered. They should furthermore respect all exact 
symmetries of the lattice model. The coefficients of these terms then have to be 
adjusted such that BRST invariance is restored in the continuum limit. A lattice 
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discretization of a nonlinear gauge, J2fi{dfiAf^ + Aj^} = 0, has been introduced in 



ref. |TI|, while a lattice discretization of the Lorentz gauge Y^nd^Afj, = was later 



given in ref. [0. It was pointed out in ref. |TT| that the lattice action of the gauge- 



fixing approach can be rewritten with Higgs fields and that on the trivial orbit the 
gauge-fixing term reduces to a higher- derivative scalar field theory. Higher-derivative 
scalar field theories have been studied recently in a series of papers [14|, but in a 
very different context. 

The Higgs field represents the longitudinal modes of the gauge field. The central 
question is then whether the fluctuations of the longitudinal modes alias the Higgs 
field are sufficiently reduced by the gauge- fixing term that the chiral nature of the 
fermions does not get spoiled |0. In order to study this question, it is useful to 
introduce the notion of a reduced model. The reduced model is obtained by keeping 
only the Higgs fields, setting the gauge field equal to zero. (In the model without 
explicit Higgs field, this corresponds to restricting the gauge field to the trivial orbit.) 
This reduced model should then have a continuum limit with free fermions in the 
desired chiral representations of the gauge group. 

In this paper, we will consider the reduced version of the purely bosonic model 
with compact U(l) symmetry that apart from the usual plaquette term includes the 
Lorentz gauge-fixing term of ref. |TB[ and a mass counterterm for the gauge field. We 
will derive the phase diagram of the reduced model, using mean-field and numeri- 
cal techniques. We will demonstrate the existence of a continuous phase transition 
between a phase with broken symmetry (FM phase) and a so-called helicoidal ferro- 
magnetic (FMD) phase, where a continuum limit can be defined. The nature of this 
continuum limit will be investigated using both perturbation theory and numerical 
simulations. In the framework of the Smit-Swift model we will demonstrate in two 
forthcoming publications |TB|, |TB[ that the unwanted Higgs-fermion bound states do 
indeed not emerge, when fermions are coupled to this model. The spectrum of this 
fermion-Higgs model contains only a left-handed fermion which couples to the gauge 
field, if we would turn it on again, and a right-handed free fermion which decouples 
from the gauge field in the continuum limit. 

The outline of the paper is as follows: In Sect. ^ we start from the full action 
of the gauge-fixing (Rome) approach and reduce it step by step to the "reduced 
model" which will be the subject of this paper. In Sect. ^ we determine the phase 
diagram of the reduced model in the mean-field approximation. The magnetization 
is calculated in the weak coupling expansion in Sect. ^. We demonstrate that the 
magnetization vanishes when approaching the FM-FMD phase transition from both 
the FM and FMD phases. Sect. ^ deals with the results of the numerical simulation. 
The numerical results for the phase diagram are presented and compared with the 
mean-field results in Subsect. |5.1|. The perturbative results for the magnetization are 



compared with numerical data of high precision in Subsect. p.2|. A brief summary 



W. Bock, M. Golterman and Y. Shamir 



3 



of our results and an outlook is given in Sect. |. 

2 The Model 

In the continuum, gauge fixing is needed in order to make the integration over gauge 
orbits well-defined. The continuum gauge-fixed action for a chiral gauge theory can 
be written as 

= s,,g{a;) + s,,p(a;; + s,,,,xa;) + s,,,^ua;;c,c) . (2.1) 

Here 

is the fermionic part of the action. Only the left-handed component of the fermion 
couples to the gauge field. S'c,g(^)!1) gauge action, S^g.fX^^) the gauge-fixing 

action and 5'c,ghost(^^; c, c) is the Fadeev- Popov ghost action. For the Lorentz gauge. 



(2.3) 



ScghosM"^] C, C) = ^ Ca [Sa,b + 9 fabc d^] Cb , (2.4) 

a,b 

where c designates the complex ghost field, ^ is the gauge-fixing parameter, g the 
gauge coupling, the gauge field and fabc are the structure constants of the gauge 
group. The continuum path integral is invariant under BRST transformations which 
replace the local gauge invariance. 

The action is transcribed to the lattice using the compact lattice link variables 

Uf,x = exp{iagA^^) e G , (2.5) 

where a denotes the lattice spacing. In the following we will set a equal to one. The 
compact link variables are elements of the gauge group G and are assigned to the 
lattice links {x,x + /t). The action on the lattice can then be written in the form, 

S = Sg{U) + Sp{U; iPl, ^r) + S^s{U) + 5ghost(t/; c, c) + S,.,.{U; V'l, ^a; c, c) . (2.6) 

Only the plaquette action Sg{U) is manifestly gauge invariant. To transcribe the 
fermion action ( |2.2| ) to the lattice one has to choose a particular lattice fermion 
formulation (like Wilson, domain wall or staggered fermions). All known lattice 
fermion formulations are in conflict with local chiral gauge invariance and, as a 
consequence, S-piU^ipL^i^R) is not invariant under BRST symmetry. This means 
that counterterms have to be added to the action so that this symmetry is restored 
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in the continuum limit. We have to consider all terms with dimension smaller than 
or equal to four, which respect all exact symmetries of the lattice model. 

In this paper we will not consider fermions and focus only on the discretization 
of the bosonic part of the action ( |2.1|) . It turns out that also the gauge-fixing part 
of the action, i.e. the combination S'g.f.(f/) + Sghost{U;c,c), should be formulated 



on the lattice such that BRST symmetry is broken. The reason is a theorem [|17 
which states that the partition function itself, as well as expectation values of gauge 
invariant observables, vanish in a lattice model with exact BRST invariance, due to 
the existence of lattice Gribov copies (see also ref. ||T8[| ). 

As a second simplification we choose U(l) as gauge group. This choice makes 
both the analytical and numerical calculations considerably easier since the ghost 
action in eq. ( |2.4| ) does not depend on the gauge potential and there are also no 
counterterms which couple the ghosts to other fields in the action. This implies 
that in the abelian case the ghost sector can be dropped completely from the path 
integral. 

Finally, as a third simplification we include only the gauge-boson mass countert- 
erm for the gauge field (this is the only counterterm of dimension two) and ignore 
all counterterms of higher dimension. The coefficient of this mass counterterm has 
to be tuned such that the photons are massless. For all dimension-four counterterms 
without derivatives it has been argued in ref. [|13| that they do not alter the phase 
structure of the physically relevant region of the phase diagram (the existence of a 
continuous FM-FMD phase transition). We believe that this remains also true if all 
other counterterms are included in the action. 

The U(l) model we are studying in this paper is then defined by the path integral, 

Z = J DU exp(-5(f/)) (2.7) 

where the action is given by 

S{U) = Sg{U) + Sg.f.(f/) + S^{U) , (2.8) 

Sg{U) = 4 E{l-Re?7^..} , (2.9) 



9 XflU 



S'g.f.(f/) = K 



\x,y,z X ) 

SUU) = -^T.{u,x + Ul} . (2.11) 

fix 

Here ?7^^^ = Ufj,xU^x+iiUj^^j^pUl^ is the usual plaquette variable, 

^{UU = E {U,.x 5x+f.,y + 5x-(.,y - 2 4,,} (2.12) 
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is the covariant lattice laplacian, 
with 

V..^ = 

2i 

and 



V, 



llX — fl 



V, 



fix 



g A^, + 0{{g A 



^X I 



(2.13) 



(2.14) 



(2.15) 



The reader can easily verify that the gauge- fixing action ( |2.10 ) reduces in the clas- 
sical continuum limit to eq. ( |2.3| ). The gauge- fixing term ( |2.3| ) can be transcribed 
to the lattice in many different ways. The choice in eq. ( p.lOp is motivated by the 



following important properties ||13|| : 

1. The action ( |2.10| ) has a unique absolute minimum at U^^ 
weak coupling expansion. 



1, validating the 



2. 



The action (p.lO|) is not BRST invariant. This is related to the fact that it 
cannot be written as a square of (a discretized version of) the gauge- fixing 



functional, I]^9^v4^. The theorem of ref. |Tj] therefore does not apply in our 
case. 



3. The action ( 2.1CI|) leads to critical behavior in the continuum limit g ^ {). 



In connection with item 1 we note that the naive discretization of the gauge-fixing 
action 



:2.i6) 



where eq. (|2.14|) was used to transcribe A^^ in ( ^73|) in terms of f/^^., does not have 
a unique minimum. It was demonstrated in ref. |1T[ that this action gives rise to a 
dense set of lattice Gribov copies. Such a dense set of Gribov copies may still give 
rise to strong fiuctuations of the longitudinal gauge degrees of freedom, a situation 
which we want to avoid from the start. 

Some first information about the phase structure of the model is obtained in the 
constant field approximation. In this approximation the lattice link field in eq. (|2.5|) 
is replaced by a constant gauge field that is independent of x. All terms which 
contain derivatives of the gauge field vanish when we insert this constant gauge field 
into eq. (|2.8|) and therefore we obtain an expression for the classical potential. After 
expanding the resulting expression for the classical potential in powers of g we find 




(2.17) 
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Figure 1: The (k, k) phase diagram of the reduced model ( ^.24 ) contains four 
different phases: a ferromagnetic (FM), an antiferromagnetic (AM), a paramagnetic 
(PM) and a directional ferromagnetic (FMD) phase. The phase boundaries in the 
mean-field approximation, obtained in the infinite volume limit, are represented by 
the solid lines. The numerical results for the phase transitions, obtained by scanning 
the parameter space in k direction on the 4'^, 6^ and 8^ lattices are marked by the 
triangles, crosses and circles. The error bars are omitted in all cases since they are 
smaller than the symbol size. 
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where the . . . represent terms of higher order in g'^. The coefficient of the term 
quadratic in has to be tuned such that the gauge-boson mass vanishes. The 
value of K where the gauge-boson mass vanishes at a given value of defines a 
critical point /tFM-FMD(fl'^)- Using eq. (|2.15|) we can also replace g^ everywhere in 



the series by k,. At tree-level we find, from eq. (|2.17| ), 



KFM-FMD(ft:) = , (2.18) 

which is a good approximation only in the large k region of the phase diagram. A 
first glance at the phase diagram in fig. |l] (which was obtained on the trivial orbit, 
i.e. Ufj,x = gxgl+fj) shows that the mean field and numerical results for the FM-FMD 
phase transition at large k are indeed very close to the k = axis. We will show 
later that the critical coupling kfm-fmd('5) is shifted by perturbative corrections to 
a small positive value. 



The minimization of eq. (p.l7|) shows that 



{g Af,) = , for k > kfm-fmd , 

{g Afi) = ± - ft^FM-FMol/ (6 k))'^^'^ , for K < kfm-fmd , 



for all /i = 1, . . . , 4 H^. This implies that k = kfm-fmd = corresponds to a phase 



transition between a ferromagnetic (FM) phase, where (^^) vanishes and the gauge 
boson has a nonzero mass, and a so-called FMD phase with a nonvanishing vector 
condensate (A^). It will become clear shortly why these phases are called ferromag- 
netic. The abbreviation FMD stands for ferromagnetic directional to express the fact 
that the vectorial vacuum expectation value (A^) induces a space-time direction and 
that, as a consequence, hypercubic rotation invariance is broken. The continuum 
limit of the model corresponds to the continuous phase transition between the FM 
phase and the FMD phase. 

To investigate the properties of the model ( p.8|) beyond the constant field approx- 
imation we can study fluctuations around the classical ground state in the FM and 
FMD phase by expanding the observables in powers of g"^, or alternatively powers 
of cf. eq. {^J^. 

Next we demonstrate that the model defined by the path integral ( p.7| ) is equiv- 
alent to a gauge-Higgs model which is manifestly gauge invariant. The lattice path 
integral in (|2.7| ) can be rewritten as 

Z = JdU exp{-S{U^x)) 

= JdU expi-SigxU,x9U)) 

= j DUDct) eM-S{ct>\U,A.+i.)) (2.20) 
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where in the second hne we have performed a local gauge transformation, U^x 
dxUnxgl+fi, Ox and made use of the gauge invariance of the Haar measure. The 
gx^s drop out of the plaquette action because it is manifestly gauge invariant. The 
third line is obtained after integrating both sides of the equation ( p.20|) over the 
gauge degrees of freedom with gx = 0^. and using the fact that / Dcj) = 1. This 
simple transformation shows that the longitudinal gauge degrees of freedom turn 
into group- valued Higgs fields. The new action S{(f)l.U^x(px+ii) is now invariant under 
the gauge transformations 

Uf,x hxU^xhl+f, , (px hx(t)x ■ (2.21) 

We will refer in the following to Sy = S{U^x) in eq. (|2.8| ) as the action in the vector 
picture and to Sh = S{(f)l.U^x(px+fi) in eq. ( |2.2U| ) as the action in the Higgs picture. 
The two actions are related by 

SviU,x) = Sh{U,x; (2.22) 
and all observables in the vector picture are mapped onto corresponding observables 



in the Higgs picture (see also ref. |jTT 



In this paper we will study a reduced model defined by the action ( p.8|) on the 
trivial orbit, Ufj^x = 9x^ dl+fi- In the Higgs picture, cf. eq. (|2.20| ), the reduced model 
is obtained by setting f/^j, = 1. The reduced model is then defined by the following 
lattice path integral 



Z = J D(f) exp(-S(0)) (2.23) 
S{<P) = j:<Pii°<P)^ + ~^T.{<Pii°'<P)^-Bl} , (2.24) 



where Bx is given by eq. ( p^.l3 ) with 



y^^- = ^ (4 0=.'+A - 0l+A ■ (2-25) 



Eq. (|2.24|) defines a higher-derivative scalar field theory. 

As a first step we will investigate in the following the phase diagram of the 
reduced model. Eq. (|2.24| ) shows that the partition function is invariant under the 
symmetry 

K, — > —K — 32 K , K, ^ K, , 4>x ^ 4>x , (2.26) 

where 

e. = (-l)^(^), E{x)=Y.x,. (2.27) 
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This implies that the phase diagram is symmetric under reflection with respect to 
the line 

k + 16k = 0. (2.28) 

For K = we recover the XY model in four dimensions whose phase diagram 
consists of three different phases: a broken or ferromagnetic (FM) phase at k > 
kfm-pm > 0, a symmetric or paramagnetic (PM) phase at kpm-am < k, < kfm-pm 
and an antiferromagnetic (AM) phase at /t < kpm-am < 0. The symmetry ( |2.26| ) 
implies that, kpm-am = — ^fm-pm- Numerically it has been found that kfm-pm ~ 
0.15. The order parameters which allow us to distinguish between these phases are 
the magnetization 

v = \{(l).)\ (2.29) 

and the staggered magnetization 

vam = \ {ex 4>x)\ ■ (2.30) 

Both quantities are not invariant under the global U(l) symmetry, and we have 
taken the modulus to eliminate the ambiguity due to the constant field mode. The 
FM phase is characterized by f > 0, vam = 0, whereas in the AM phase v = 0, 
Vam > 0. Both order parameters vanish in the intermediate PM phase. 

As explained above, at large R we expect to find a new phase transition between 
the FM and the FMD phase, which at tree level is given by eq. (|2.18 ). (In the 



following we will retain the name FMD also for the reduced-model version of the 
FMD phase.) The FMD phase is characterized by a new vector order parameter g^, 
< < 271, which is nonzero in the FMD phase and vanishes in the FM phase. It 
is equal to (vr, vr, vr, vr) in the AM phase. As a generalization of v and vam we define 
a helicoidal magnetization 



X 



(2.31) 



which is nonzero in the FMD phase. It is easy to see that reduces to v in the 
FM and to wam in AM phase. 

When ignoring fluctuations around the ground state, the vector field V^^^ cf. 
eqs. ( |2.14| ) and ( p.25| ), in the FMD phase is given by 



V,x = vi:,q, + 0{([1) , (2.32) 

showing that plays the role of the vector condensate in the reduced model. We 
mention in passing that phases with nonvanishing have been intensively investi- 
gated in lower dimensions in condensed matter physics and are known as helicoidal- 
ferromagnetic phases (see ref. [|T^ for a recent review). 
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To further substantiate the statements about the phase diagram made in this 
section we will determine in the next section the phase diagram of the reduced model 
( |2.24| ) in the mean-field approximation. Numerical data for the phase diagram are 
presented in Sect. |5.1| and compared with the mean-field results. 



3 The Phase Diagram in the Mean-Field Approx- 
imation 

In the following we will perform a mean-field analysis of the phase diagram in d 
dimensions. 

A central problem of the mean-field approximation in more complicated ferro- 
magnetic systems is the choice of the mean-field ansatz which in a given region of 
the parameter space leads to the absolute minimum of the free energy. Usually there 
exist many different choices and it is not straightforward to pick an ansatz which 
leads to the absolute minimum of the free energy. Based on the discussion of the 
previous section we decided to consider the ansatz 



f exp Xi}j , (3.1) 



where Qfi, < qfj_ < 2-71 are real phases and if plays the role of a magnetization. 
Depending on the value of this ansatz can distinguish between phases with fer- 
romagnetic {q^ = 0, fi = 1, . . . ,d), {p = V, antiferromagnetic {q^ = n, fi = 1, . . . ,d) 
ordering, (p = f am, and phases with a helicoidal magnetization (g^ ^ 0, vr, for at 
least one component fi), = v^. Similarly, we take for the magnetic field the 
ansatz 



= h exp ' (3-2) 



where h is the mean-field magnetic field. 

Using eqs. ( |3.1|) and ( p.2|) and following the steps of the standard mean-field 
calculation (see for example ref. pOf), we obtain for the free energy of the reduced 
model ( p:2lD 



(3.3) 



where L is the extent of the lattice in spatial and temporal directions, 

/(i)(?; K,K) = -2{4dK + k) F{q) + 2 ^ {2F{qf - d) + ^ F{2q) {2d + 1) , (3.4) 

Id 
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f(^)(q-R,^) = -^ {6F{2qf-4{3d+l)F{2q) + d{10d-l)) , (3.5) 
R, k) = (2 Fi2qf - 2 (d - 1) F{2q) - d) , (3.6) 

/(^)(g;^,«:) = -^(F(2g)-d)^ (3.7) 



and 



^(9)=Z!cosg^ (3.^ 



1 

IqQi) = — da exp(±/i cos a) (3.9) 

TT JO 



is the modified Bessel function of zeroth order. We have dropped in eq. ( |3.3| ) all 
terms that depend neither on nor on q. The saddle-point equations read 

^ = L''S^2h + 2 I f^\q- ^, «:)| = , (3.10) 

L"^' sing, |^y,2-2^?«(g;/?,«:)| = 0, (3.11) 



'dh /g(2/i) 



-2L^L-hm\=0, (3.12) 



where 



gl!\q;K,K) =2{AdK + K)-8K F{q) - ^ (2 d+ 1) cosg, (3.13) 



gj^\q;li,^) = --{12 F{2q)- 12 d~ A) cosg, (3.14) 



K 



,(3) 



16 



gl!\q;K,K) = KiFi2q)-^^id-l)) cosg, (3.15) 
gl'^Xq-.K^K) = I (F(2g) -rf) cosg, , (3.16) 



and 



/.(A) = (3.17) 

is the modified Bessel function of first order. 

From these d + 2 equations we can compute the d + 2 fields (p, q^j, and h as 
functions of the parameters H and k. The phase boundaries are defined as the lines 
in the (k, k) parameter space where various combinations of the order parameters 
(f and q^ vanish. 
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The variable h can be eliminated from the saddle-point equations in regions where 
ip is very small (which is the case close to the PM phase, where ip vanishes) and the 
ratio /i(2/i)//o(2/i) in eq. ( |3.12|) can be expanded in powers of h, Ii{2h) / Io{2h) = 
h + 0{h^). 

Usually there does not exist a unique solution of the saddle-point equations in 
a certain region of the parameter space. It is therefore important to substitute the 
various solutions back into the expression for the free energy ( p. 3D and to pick out the 
solution that corresponds to the absolute minimum. In practice it can happen that 
certain phases remain undetected because the mean-field ansatz was too simple. In 
the following we will consider also another ansatz to search for a ferrimagnetic (FI) 
phase in a certain region of the parameter space which cannot be probed with the 
ansatz ( |3.1| ). Because of this uncertainty of the mean-field calculation it is important 
to determine the phase diagram also numerically. 

We furthermore note that the free energy in eq. p.3| ) is invariant under the 
symmetry ( p.26| ), 

— s> TT — , K —K — 8 dR , k — k , (3.18) 

which implies that also the phase diagram in the mean-field approximation is sym- 
metric (but for the interchange vr — g^, which maps the FM onto the AM 
phase, etc.) with respect to the line k + AdH = which in four dimensions turns 
into eq. (^^. 

In the following paragraph we present our mean-field results for the phase bound- 
aries and briefly explain how they were obtained: 

• FM-PM and PM-AM transitions: The transition between the FM and 
PM (PM and AM) phases is obtained by approaching the transitions from 
within the FM (AM) phase where F(g) = d, F{2q) = d (F(g) = -d, F{2q) = 
d) and p = v {ip = wam) approaches zero. The relation for the FM-PM (PM- 
AM) phase boundary is obtained by expanding eq. (|3.10| ) in powers of p and 
equating the part that is linear in p with zero. Eq. ( p.ll[ ) is trivially fulfilled 
because sin g^ = in the FM (AM) phases. The FM-PM and PM-AM phase 
boundaries are respectively given by the relations, 

1 31 _ 

«;fm-pm = 2rf~*^^^~^''"'*32'^' {'i-l^) 

1 /194 , 31\ , , 

«:PM-AM = -^-(^^-^j (3-20) 

The reader can easily verify that the two solutions are related to each other 
by the transformation ( p. 18 ). The FM-PM and PM-AM phase transition lines 

. The corresponding solution of the 



intersect aX k = Ki 



" V 16 8 " 
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saddle-point equations minimizes the free energy ( ^.3| ) only in the interval 

16 

[d{2d+ 1)]-' = 1^2 , (3.21) 



and the two straight lines ( |3.19|) and ( p.2(]| ) therefore form the boundary of 



the PM phase only in that interval. 

FMD-PM transition: The magnetization = f h vanishes when we ap- 
proach the phase boundary from the FMD side which means that eqs. ( |3.10| ) 
and (p.ll|) can be expanded in powers of (f. Unlike in the previous case 
sin Qfj, ^ for at least one component fi = 1, . . . ,d and hence —d < F{q) < +d. 
From the term in eq. (|3.1CI|) that is proportional to if we obtain 



2(4dK + K) Fiq) -2 k (2 Fiqf - d) F{2q) (2 + 1) = 1 . (3.22) 

16 



After summing the i = 1 term inside the curly brackets in eq. (|3.11 ) over 
we obtain for F[q) the formula 

, A dn + K , ^ 

and after first multiplying the same term with cosg^ and then summing over 
/i we obtain an expression for F{2q), 

^{2d+l) F{2q) = {AdK + K) F{q) -4 k F{qf - ^ {2 d + 1) d . (3.24) 
16 16 

After substituting these two solutions into eq. ( p.22| ) and after a few trivial 
algebraic manipulations we obtain the following solution for the FMD-PM 
phase transition: 



^FMD-PM — ~4 d K ± 



Wd 



l- — d{2d + 33) 
16 



(34rf+l). (3.25) 



This solution describes an ellipse located around the symmetry axis, k+4 d k 
0. Eq. ( p.25| ) forms the boundary of the PM phase in the interval 



< K < 16 [rf (2 rf + 33)]" = K3 , (3.26) 

as the corresponding mean-field solution does not lead to an absolute minimum 
of the free energy ( p. 3D in the region where k <K2- 
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FM-AM phase transition: Above we pointed out that the FM-PM and 
PM-AM phase transitions intersect ai k = Ki. This means that the FM and 
AM meet at this value of k. Two different scenarios are imaginable for the 
phase structure in the region k < ki: 1) The FM and AM phases meet at the 
4 d K, + K, = symmetry line with the magnetizations v and f am exhibiting a 
jump at this line, or 2) the FM and AM phases are separated by a FI phase 
in which both order parameters v and f am are simultaneously nonzero. The 
mean-field ansatz ( |3.1|) is not suited to detect such an intermediate FI phase 
since cannot be equal to and tt at the same time. We therefore calculated 
the free energy also for the ansatz of the form 

(px = v + vam ex , (3.27) 

which allows us to probe for a FI phase. Our calculation however shows that 
scenario 2) leads to a larger value of the free energy. Also our numerical data 
in four dimensions give clear evidence for the correctness of the first scenario. 

FM-FMD and FMD-AM transitions: The FM-FMD (FMD-AM) phase 
transition is characterized by — > (g^ it) for all /i = 1, . . . , c?. It is 
difficult to determine the location of these phase transitions analytically be- 
cause ip does not vanish at these two phase transitions and eqs. ( |3.1CI| )- (|3TT^) 
cannot be expanded in (p. We will show in the next section that the magneti- 
zation actually vanishes at the FM-FMD (FMD-AM) phase transition. This 
phenomenon is connected to the infra-red behavior of the higher-derivative 
action and therefore cannot be understood in the framework of the mean-field 
approximation. 

The FM-FMD phase transition is determined by eliminating the fields h and 
ifi from equation (|3.12|) and 



-4:dK-4:dR{2d+l) + -d(2d+l) + — Kd (-4(1 + 5)^'^ 

8 16 

--Kdip^ + 2- = 0, (3.28) 

8 if 

2 K - ^ (2 + 1) - i K ^2 ^ ^ (rf + 1) (/?^ = . (3.29) 



Eqs. ( |3.28|) and ( p.29| ) are obtained from eqs. ( |3.10| ) and ( |3.11| ) in the limit 



g — * 0. Note that eqs. (|3.12|) , (|3.28| ) and (|3.29| ) can only hold simultaneously 



on the FM-FMD transition curve, where F{q) = F(2q) = d. The location of 
the FM-FMD transition can be calculated analytically in two special cases: 
at the PM-phase boundary, the terms in eq. ( |3.29| ) which are quadratic and 
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quartic in (p can be ignored, and the FM-FMD transition is given by the 
intersection of the PM-phase boundary and the hne 

«;fm-fmd = g (2 + 1) . (3.30) 

Similarly the FMD-AM phase transition is given by the intersection of the 
PM-phase boundary and the line 

kfmd-am = -8 d K - - {2 d + 1) . (3.31) 

o 

The FM-FMD phase transition can also be calculated analytically in the limit 
oo. Eq. ( lOSi ) implies that h = Rd (128 d + 63)/32 + 0(1) for k ^ oo, 



and the expansion of (p, cf. eq. ( p.l2|) , gives (p = 1 — l/(4/i) + 0{l/h ). After 



substituting these two formulas into eq. ( |3.29| ) we find 

2 (4 + 3) 

«:fm-fmd - ^^28 d + Q3) d ' " ^ ' ^^'^^^ 

In four dimensions we have determined kfm-fmd(j^) also at a series of inter- 
mediate K values by solving eqs. ( p.l2|) , ( |3.28|) and (|3.29| ) numerically. The 



results for fi:FM-FMD('^) are listed (with an accuracy of four decimal places) in 
the second column of table |l|. We have also included the mean-field value of 
evaluated at k = /tFM-FMD('5) (third column). The corresponding numerical 
values of the FMD-AM phase transition can be easily obtained from the data 
in the second column of table ^ by making use of the symmetry (|3.18|) . 

The mean-field phase diagram for d = 4 is displayed in fig. |l]. The phase bound- 
aries are represented in this plot by the solid lines. The FM-FMD phase transi- 
tion line was computed by solving eqs. ( |3.12| ), ( |3.28| ) and ( p.29| ) numerically. The 



FMD-AM transition was computed from the the FM-FMD phase transition data 
by making use of the symmetry ( p. 18] ). The various symbols in fig. |I] represent the 
results of the numerical simulation and will be explained later. 



4 Weak Coupling Expansion 

Eq. (|2.15| ) suggests that the weak coupling expansion should be performed in 1/k. 
We consider small fluctuations of the fields around the classical ground state, 
which we parametrize by the Goldstone field 6x'- 

4>x = exp ^q^,x^ + i 6^/ V2Rj . (4.1) 
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mean- 


field 




approximation 




/^FM-FMO 




0.013 


0.01467 


0.21887 


0.014 


0.01566 


0.42089 


0.015 


0.01597 


0.52455 


0.020 


0.01639 


0.73962 


0.025 


0.01621 


0.81831 


0.030 


0.01609 


0.85916 


0.050 


0.01608 


0.92369 


0.1 


0.01626 


0.96377 


1 


0.01649 


0.99651 


10 


0.01651 


0.99965 


100 


0.01651 


0.99997 


oo 


0.01652 


1.00000 



Table 1: The mean- field results for the critical coupling kfu-fmd{^) o,nd the order 
parameter evaluated at n = kfm-fmdI^) o'^e given in columns two and three for 
several values ofk. The value o/^fm-fmd at H = oo was calculated from eq. ( ^.3^ ). 
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This ansatz holds both in the FM phase, where = and also the FMD phase, 
where ^ 0,7i for at least one component = 1, . . . , 4. The phases in the FMD 
phase can be computed by minimizing the effective potential, or alternatively could 
be taken from the numerical simulation. 

To calculate v in perturbation theory, we first insert eq. ( [4.1| ) into the action 
(p.24|) of the reduced model and expand it in powers of 9^, 



s = \l e{-k) A-\k) 0{k) 



+ / v,(A;i, k2, ks, k,) e{ki) e{k2) e{k^) e{k,) + . . . , (4.2) 

J ki,k2,k3,k4 

where 

f f^^ d^k , , 

Iri 12^- <^'^) 

The subscript q indicates that the propagator Aq(fc) and four-point vertex function 
Vg{ki, k2, k^, k^) depend on the phases q^. The inverse propagator A~^{k) is given 

by 

= (Ecosg^2 (l-cosA;^)| -8 - 4) ^ cos (1 - cos A;^) 

+ < 2 E sin sin fc^ ^ + 4 X! ^^^^ sin^ q^, (1 - cos k^) 

- I'^sin 2q^ sin fc^ > ~ 2 ^ sin^ q^ cos^ q^ sin^ fc^ 

I M ) fiu 

cos 2 (1 — cos kfj) , (4.4) 

with 

= - . (4.5) 

K 

(The expression for the propagator simplifies if g is a solution of the classical saddle- 
point equations, cf. eq. ( [4.13|) below. As we will discuss below, however, this is not 
in general the case in finite volume.) At tree- level, the FM-FMD phase transition 
line is given by = 0, cf. eq. ( p.l8|) . When approaching the FM-FMD phase 
transition line, eq. (W^ ) reduces to 

Aq{k) oc l/ik^f , K kfm-fmd (4.6) 



for small k. The propagator ( [4.6|) leads to infra-red divergences. A similar situation 
is encountered in two space-time dimensions where infra-red divergences occur for 
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massless bosons with an ordinary kinetic term These infra-red divergences 



are not only an artifact of the tree-level propagator, but occur in the full theory 
when the continuum limit k — > kfm-fmd is performed. The qualitative agreement 
with the two-dimensional behavior will be demonstrated below both analytically and 
numerically. Here we note that the situation we encounter in the reduced model is 
similar to the situation of the XY model in two dimensions. The FM-FMD phase 
transition line behaves like the spin-wave phase, where critical exponents depend 
continuously on the coupling constant. Below we will show that the magnetization 
V (helicoidal magnetization vu) vanishes oc — kfm-fmdI''*-'^'' when k \ kfm-fmd 
(k y /tpM-FMD) and that the critical exponent rjiji) depends continuously on k. 

It is useful to distinguish between observables which are invariant under the 
global U(l) symmetry (symmetric), like the two-point function and others 

which are not invariant (nonsymmetric), like the magnetization {(px)- For symmetric 
observables, the weak coupling expansion is infra-red finite, because all interactions 
involve derivatives. The situation is different for nonsymmetric observables, such 
as the magnetization. The real expansion parameter is not 1/k in this case, but 
(logm^)/fi;. This means that in order to obtain a nondivergent result in the limit 

— 0, one should perform a resummation of infinitely many diagrams. 

Using eqs. ( [4.1D and ( [4.2] ) the magnetization, cf. eq. ( |2.31| ), can be calculated to 
one-loop in perturbation theory. 



[ Aq{k) + higher order corrections . (4-7) 

A K Jk 

The integral in eq. ( [4.7| ) is infra-red divergent in the limit — > and, as mentioned 
in the previous paragraph, in order to obtain a finite result we have to resum the 
higher order diagrams (with two and more lines) in fig. ^ which arise from the 
terms proportional to 6^"' in eq. ( [4.1|) with n > 1. This resummation of diagrams 
gives 



exp ( — [ Ag{k) \ + higher order corrections (4.^ 




V 



+ higher order corrections , (4-9) 

where we used that at tree-level kfm-fmd = 0. The higher order corrections in 
eqs. ( |4.8| ) and ( ^.9| ) are due to quartic and higher order interactions which we have 
ignored. The critical exponent t] in eq. ( [4.9|) is given by 



"=64^- <"°> 

Eq. ( [4.9|) shows that the magnetization vanishes at the FM-FMD phase transition 
with a critical exponent t] that depends on R and differs from the Gaussian exponent 
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(b) 



Figure 2: Feynman diagrams for the magnetization (a) and the scalar field propa- 
gator (h). 

1/2 of the XY model at /5 = 0. (In order to show that also the helicoidal magne- 
tization fH vanishes for k /^fm-fmd with the critical exponent ([4.10|) , one uses 
eq. (|4.4| ) where q is the nontrivial solution of eq. (|4.13|) below.) 

Using relation ( [4 .81 ) and replacing the integral by a lattice sum, we can compute 
f H on a finite volume at any value of H and k in the FM and FMD phases provided 
that we know the function q = q{R, k). On a finite lattice with periodic boundary 
conditions for the scalar fields the phases q^ can only take the values 



q^ = 2n^7r/ L 



n, 



0, 



L-l 



(4.11) 



where L designates the extent of the lattice in temporal and spatial directions. The 
g's can be calculated at tree level by minimizing the classical action 



So{q- K, k)=L^ { -2 (16 k + k) F{q) + 4 k F{qf 



sin^ qf. 



(4.12) 



where is the volume of the lattice and F{q) is defined in eq. (|3.8| ). We note that 
after expanding eq. ( |4.12| ) in powers of q and setting q^ = gA^, we recover to leading 



order the classical potential in eq. (|2.17| ). In infinite volume, the phases q can be 
determined from the four saddle-point equations. 



2 (16 K + k) — 8 K F{q) — A k ^ sin^ q^ \ cos q^ 



sm q^ 







(4.13) 
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In finite volume, the value of q that minimizes So{q] k, k) at a given {k, k) point in the 
FMD phase is determined by computing So{q', k, k) for all vectors q numerically, 
and picking out the ones with the smallest value of So{q] k, k). (The minimum will 
respect the lattice symmetries and hence will in general not be unique.) From the 
resulting g- values we have computed the observable F{q), cf. eq. (|3.8| ), for several 
(k, k) points in the FMD phase. The discussion of this quantity will be postponed 
to Sect. |5]^ where we will compare it with the results of the numerical simulation. 

We will from now on focus on the physics in the FM phase where g = 0, and 
calculate the magnetization to one higher order in perturbation theory. From this 
calculation we will obtain another estimate for the critical coupling kfm-fmd('5)- 

The vertex function Vo(fci, A;2, ^3, k^) in eq. ([4.2|) is given by 



Vo(fci, ^2, k^) 



1 

4 K 



^(1 - cos(fci + k2) 



^ 3 



^(1 - cosfci^) 



+ ^ sin ki^ sin k2^ sin sin k^^ + ~ cos(/ci + ^2) 

^(1 - COS ki^) \ 6{ki + k2 + k3 + k^) , 



(4.14) 



where the first three terms arise from the kX^xI^K'-'^^)^ ~ -^x) term and the two 
terms proportional w? from the —kY1,x 0i('^0)x term in the action ( p. 241 ). We note 
that the vertex function ( [4.14| ) can be easily rewritten such that it is symmetric with 
respect to the momenta fci, A;2, ^3 and ^4. For the perturbative calculation it does 
not matter which form is used. 

After carrying out the higher-order calculation, the magnetization can be written 
in the form 



,1— loop 

(fc) 



exp ( — [ Aoi_ioop(^)l + higher order corrections , (4.15) 

\ A K Jk ' J 

Ao(fc) 1 



1 + Ao(A;)S(A;) Aq '(A;) + E(A;) 



where Ao(A;) is the tree-level propagator in the FM phase (cf. eq. 
and S(fc) is the self-energy. 



(4.16) 
with g = 0) 



S(A;) 




^(1 - cos(p + A;)^ 



-\ 2 



-\ 2N 



2 5](l-cosA;^) 



+ 



cos 
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- sin^ sin^ Pi/ - 2 ( XI sin k^ sinp^ j 



(4.17) 



The propagator (^4.16 ) already involves a resummation of diagrams shown in fig. 03. 
The magnetization in eq. (|4.15| ) has been obtained by performing the resummation 
of diagrams in fig. ||a, but now using the propagator Ao,i_ioop(A;) instead of Ag(fc). 
In order to compare the perturbative formulas ( |4.8| ) and ( ^.15| ) with the results of 
the numerical simulations (see Sect. |5.2| ) we have to evaluate the lattice integrals 
in eqs. ( ^.8| ) and ( |4.15| ) on a finite lattice. The integrals are replaced by sums over 
the lattice momenta. In this context we note that these finite lattice sums do not 
include the zero mode, k = 0. The zero mode decouples from the action, and gives 
rise to a phase which disappears after taking the modulus in the definition of v and 
vn in eqs. (|2:29|) and ( ^^311) . 

The critical coupling can be calculated by expanding Ao,i-ioop(A;)^^ for small 
momenta in powers of k'^, 



Ao,i-ioop(A;)-' = Ao(A;)-^ + S(A;) = a{K, k) {k') + b{K, k) {k^ + 
and equating the coefficient a(fi;, k) to zero. 



(4.18) 



a K, K 



= ^ + 4 /|| ^(1-cosA;^) ^ cos - ^ sin^^ 

K K Jk [ ^ fl 



cos^ k,, 



Ao(A;) = . 



(4.19) 



This leads to the one-loop estimate 



fiFM-FMD('«) = - \ \ 51 cos A;^ - ^ sin^ A;^ I 

■^^ [ fl V fl ) 



X 



1 -2 



y^(l — cos k 
0.02993 , 



(4.20) 
(4.21) 



which is about a factor two larger than the mean-field value, cf. table |l|. 



5 Numerical Results 
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5.1 Phase Diagram 



To simulate the reduced model defined by the path integral ( p.23| ) we have imple- 
mented two different Monte Carlo algorithms, a five-hit Metropolis and a Hybrid 
Monte Carlo algorithm. The results for the various observables agree nicely within 
the AM, FM and PM phases. We find however that the Hybrid Monte Carlo algo- 
rithm gets much more easily stuck in metastable non-equilibrium states in the FMD 
phase. We therefore have generated the bulk of the data presented in this paper 
with a five-hit Metropolis algorithm. 

To map the phase diagram we have measured the following observables: 



The magnetization 



(5.1) 



which is the order parameter for ferromagnetism and 
the staggered magnetization 

" 1 



(5.2) 



which is the order parameter for anti-ferromagnetism. 

The helicoidal magnetization 



exp 



was used to map the FMD phase, where the four real phases g^, /z = 1, 
were determined for each configuration from 



= Im Log 



(5.3) 
...,4 

(5.4) 



Apart from these quantities we have also measured the internal energy density 

o / 1 



4L4 



-'x+fi 



(5.5) 



fix 



of the mass counterterm and 
the quantity 

c(g/.) = (cosg^) , (5.6) 
where the phases were calculated for each configuration by means eq. 
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Figure 3: Scans in k direction on a 4^ lattice at k — 0.2 (fig. a) k — 0.1 (fig. h) 
K — 0.02 (fig. c) and k — —0.03 (fig. d). The results for v, v^m o-nd v^i are 
represented in the plots by the circles, triangles and crosses. The error bars are 
omitted because they are in all cases smaller than the symbol size. The dashed lines 
in figs, c and d are obtained in the mean-field approximation. 
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We have taken the modulus in eqs. ( p.l| - p.3| ) for each configuration, because, in a 
finite volume, the constant field mode gives rise to a slow drift of the magnetization 
through the group space. (Taking the absolute value is a standard method which 
allows us to avoid the introduction of an external magnetic field.) 

The discreteness of the phases ( [4.1f| ) poses a problem for the simulation in the 



FMD phase because each transition of one q to another q behaves very much like 
a first order phase transition and hence is accompanied by metastabilities. We find 
that these metastabilities become more severe when the lattice size is increased. 

To determine the phase diagram we kept the parameter k fixed and performed 
simulations at a large number of n values. Each of these vertical scans, cf. fig. 
has been started in the FM phase. We lowered k in fixed steps and used the last 
configuration of a run as the initial configuration at the next smaller value of k. At 
each point we skipped 10^ sweeps for equilibration and performed 10^ measurement 
sweeps. The error of an observable O was determined using the relation 



A(0) = A(0)stV2^, (5.7) 

where A(0)st is the standard deviation and Tint designates the integrated autocor- 
relation time, defined as Tint = Sa* r(At)/r(0) (see for example ||22|). The quantity 
r(At) = (O(t)O(t + At)) - (0)2 is the autocorrelation function. 

In figs, ^-d we have displayed the numerical results for the order parameters v 
(circles), vam (triangles) and vu (crosses) for four exemplary scans on a 4"^ lattice. 
Fig. ^a shows that there are three different phases at /5 = 0.2, an FM phase at 
— 0.43, where f = fn > 0, vam = 0, an AM phase at — 5.97, where 
■^AM = Vu > 0, V = and the FMD phase at intermediate n in which fn > 0, 
V = Vam = 0. The helicoidal magnetization fn exhibits jumps at certain k values 
within the FMD phase. These jumps in vn occur because the phases q^ can change 
only in discrete steps and hence have to be considered as a finite volume artifact. 
In Fig. ^ we have plotted the quantity c(g^), = 1, ... ,4 for the same scan at 
K = 0.2 as a function of n. A comparison of figs. ^ and ^a shows that the jumps in 
fH occur at the same n values where one of the components c(g^) exhibits a jump. 
All g^'s are zero in the FM phase. At k, ^ —0.43 the first component of q condenses 
(dotted line) and becomes equal to 27i/L = 7i/2. The next jump occurs when also 
the second component of q becomes equal to 7r/2 (solid line). The k value where 
finally all values of q are equal to 7r/2 coincides nicely with the symmetry point, 
K = —16k = —3.2, cf. eq. ( [3.181 ). We note that the order in which the jumps occur 
is arbitrary (because of hypercubic symmetry). The jumps at k < — 16/t follow a 
similar pattern, with the g^'s jumping from tt/2 to vr. We will show in Sect. ^.2| that 
the complicated ^-dependence of vh in the FMD phase can at least qualitatively 
be explained by the one- loop formula (|4.8|) . In fig. we have plotted the internal 
energy density function of k. In the mean-field theory this quantity is given 
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Figure 4: The observahles c^q^) (fig. a) and (fig. h) as a function ofuforR = 0.2. 
The lattice size is 4^. The four observables c^q^), 11 = 1, . . . , 4 are distinguished by 
the line type (solid, dots, long and short dashes). The error bars are omitted in both 
figures because they are not larger than the symbol size in fig. b and also not much 
bigger than the line width in fig. a. 
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by v^j J2fiCosq^ and since = 0(1), we expect this quantity to jump whenever q 
changes. The comparison of figs, ^a and ^ shows that this is indeed the case. 

A different situation is encountered in fig. ^ which shows the result of the scan 
at K = 0.1. Besides the FM, AM and FMD phases we find now clear evidence 
for a PM phase (where v = vam = vn = 0) which, as predicted by the mean-field 
calculation, extends into the FMD phase. The graph shows that four different phase 
boundaries are crossed when k is lowered from the FM to the AM phase. The small 
peak at K = —16k ~ —1.6 appears to be a finite size effect because it becomes 
smaller when the lattice size is enlarged. 

The FMD phase gradually disappears when k is lowered further. The situation 
at K = 0.02 is depicted in fig. The FMD phase has now completely disappeared, 
and the only three phases we are left with are the FM, PM and AM phases. We find 
that the PM phase extends down ton^ —0.02. The result of the scan at k = —0.03 
is displayed in fig. ^d. It shows that the phase transition between the FM and AM 
phase coincides with the symmetry line n = —16k and is obviously of first order. 
Both the internal energy density and the order parameters v and f am exhibit a 
gap at this phase transition. This gap grows from zero to one when one follows the 
symmetry line k = —16k from the triple point where the FM, PM and AM phases 
meet to k = — oo. The two dashed lines in figs. ^ and d represent the mean-field 
result for the magnetization = v in the FM and for the staggered magnetization 
V = "^AM in the AM phase, which we obtained by solving the mean-field equations 
( p.lO| ) and ( p.l2| ) for q = (0, 0, 0, 0) and q = (vr, vr, vr, vr) numerically. 

We have read off the positions of the various phase transitions from plots like 
the ones depicted in figs, ^-d and then compiled them in the k-k phase diagram 
plot in fig. |I|. The triangles were obtained on a 4^ lattice, the crosses represent the 
phase boundaries on a 6^ lattice and the circles mark the phase transitions points 
on an 8^ lattice. The solid curves are the mean-field results which we discussed in 
Sect. 1^. The plot shows that the numerical estimates for the FM-AM, FM-PM, PM- 
AM, FM-FMD and FMD-AM phase boundaries agree nicely with the mean-field 
prediction. The agreement seems to be worse for the FMD-PM phase transition. 
The numerical data indicate however that both the horizontal and vertical width 
of the PM phase shrink when the lattice size is increased and that the numerical 
results could come out closer to the mean-field result for larger volumes. 

In fig. ^we have displayed the magnetization w as a function of k for k = 0.2 and 
a series of different lattice sizes. Again, we have lowered k in small steps, and used 
the last configuration of a run as initial configuration for the next run, skipping 10^ 
sweeps for equilibration. 

First we discuss the results on the smaller lattices of size 3^ (filled circle), 4^ 
(filled square) and 5^ (filled triangles). The magnetization v exhibits a jump on 
these lattices. The helicoidal magnetization uh is identical with v in the FM phase 
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Figure 5: The magnetization v as a function of k, atK — 0.2 on 3^ (filled circle), 
(filled square), 5^ (filled triangle), 6^ (crosses), 8^ (open circle), 10^ (open triangle) 
and 12"^ (open square) lattices. 
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and (unlike v) remains of 0(1) when crossing the FM-FMD phase transition towards 
the FMD phase (the data for wh are not included in fig. It can be seen that the 
curves for v bend in more strongly when the lattice size is increased, indicating that 
V scales to zero at the FM-FMD phase transition as predicted by the weak coupling 
expansion in Sect. We will demonstrate in the next subsection that the data for v 
in the FM phase are nicely consistent with the perturbative formula ( 4.15| ) according 
to which f \ in the limit k \ kfm-fmd and L — > oo. The plot also shows that 
^FM-FMD increases with increasing lattice size. 

On the larger lattices we encounter a different behavior. The magnetization first 
starts to bend over when k is lowered but then instead of jumping to the FMD phase 
continues to decrease slowly. The jump to the FMD occurs finally at a large negative 
value of K. The jump on the 6^ lattice (crosses) for instance occurs at k, ^ —0.9 
and not at k ^ —0.2. The reason for this effect could be related to the fact that 
the FM-FMD phase transition behaves very much like a first order phase transition 
as the magnetization drops to zero very rapidly at kfm-fmd- On the larger lattices 
the fluctuations in the internal energy are very small and hence the probability for a 
jump across a large energy barrier becomes strongly suppressed. This situation did 
not change after increasing the statistics by one order of magnitude which means 
that the probability for a transition to occur at larger k is very small. In this 
context we also note that on larger lattices the system in the FMD phase ends up in 
different states when using different starting configurations and transitions to other 
states occur very rarely or not at all. The results for the various observables in the 
FMD phase are independent of the initial configuration only on the smaller lattices 
of size 3"^, 4^ and 5^. (In all other phases, our results are independent of the initial 
configurations on all volumes.) 

Fig. ^ shows that the region in the FM phase where the magnetization starts 
to bend over is shifted in all cases towards larger values of k, when the lattice size 
is enlarged. We will show in the next subsection that this finite size behavior is in 
nice agreement with the perturbative formula given in Sect. ^. This good agree- 
ment between the numerical data and lattice perturbation theory lead us to identify 
^FM-FMD on the larger lattices, i.e. for L > 6, with the point where the slope in v is 
largest, and not with the point at large negative k, where v exhibits the jump and 
fH becomes different from v. All phase transition points on the 6^ and 8"^ lattice 
which are included in fig. |l| were obtained with this criterion. 



5.2 Comparison with Perturbation Theory 

In this subsection we compare the simulation results for v in the FM and in the 
FMD phase with the perturbative formulas which we derived in Sect. ^. In fig. ^we 
have plotted once more the v and vu data of fig. ^a which were obtained at k = 0.2 



W. Bock, M. Golterman and Y. Shamir 



29 




Figure 6: The magnetization v (squares) and the helicoidal magnetization fn 
(crosses) in the FM and FMD phases as a function of k for k = 0.2. The same 
data were presented already in fig. \^a. The circles were obtained by numerically 
evaluating the perturbative expression (^~3^). The phases = q^{n, k) were obtained 
by minimizing eq. U-^W - 
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on a 4^ lattice. The circles in fig. ^ were obtained by evaluating the one-loop formula 
for f H ( |4.(j| ) numerically on both sides of the FM-FMD phase transition on the same 
lattice and at the same values of k where we performed the numerical simulations. 
The phases k) at a given value of k, were determined analytically by minimizing 
the classical action in eq. ([4.12| ) with L = 4. Fig. ^ shows that the numerical 
results are nicely reproduced by the analytic formula (4^) at large values of k. At 
smaller k values the deviations start to become larger. While we do not understand 
this phenomenon in detail, we believe that it may be related to the metastabilities 
mentioned in Sect. (Higher orders in perturbation theory could also be sizable 
though.) The jumps at which the components of q condense occur slightly delayed. 
This distorts the K-dependence of the propagator, which depends also explicitly on 
K and not only implicitly through the phases k). From the minimizing phases 
q{R, n) on the 4"^ lattice we have computed the function F{q), cf. eq. (|3.8|), which we 
plotted in fig. ^ (dashed line) versus n. The numerical results for z"^, which in the 
mean- field approximation are proportional to F{q) are represented in this graph by 
the crosses. It can be seen that the gap between the two curves widens up when k is 
lowered, which is presumably due to the systematic delay of the transition events in 
the numerical simulation. We should observe a smaller shift on a smaller lattice if 
this scenario is correct. Fig. ^ shows that the gap between the two curves shrinks 
indeed on the 3^ lattice. We also checked that the numerical metastabilities on a 
5^ lattice get stronger making the agreement with the analytic results worse. Note 
that the discrepancies between perturbation theory and numerical data occur at the 
same locations in figs. | and 0a. 

In the following we will discuss only the FM phase. The numerical metastabilities 
mentioned in the previous paragraph have an effect on the simulation results only 
in the near vicinity of the FM-FMD phase transition. 

In fig. § the magnetization v for k = 0.1 is displayed as a function of k. The lattice 
size is 8^. At each k point we have accumulated a statistics of 10^ Metropolis sweeps. 
The magnetization was measured after each sweep and its error was computed by 
means of eq. ( |5.7| ). To compare the numerical data with the perturbative formulas, 
we have numerically evaluated the integrals (replaced by lattice sums) in eq. ( [4.8| ) 
and (|4.15|) for a large number of k, values in the FM phase on the same lattice which 
we used in the numerical simulations. The obtained results for v in eqs. ( |4.8| ) and 
(|4.15|) are represented by the dashed and solid curves. Fig. | shows that the two-loop 
formula ( [4. IS]) provides, as expected, a much better description of the numerical data 
than the one-loop formula (|4.8| ). The fact that perturbation theory in 1/k remains 
valid down to such small values of k is because the actual expansion parameter is 
not I/k, but l/(167r^K) where the factor l/(167r^) comes from the loop integrals. 

The K-dependence of the magnetization v is shown in fig. ^ where we have 
plotted the magnetization data for five different k values versus k. The lattice size 
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0.4 - 



0.2 - 



Figure 7: The quantity \F{q) (dashed line), obtained from the minimization of the 
classical action and the internal energy density (crosses), obtained from 

the numerical simulation, as a function of k for k = 0.2 on 4^ (fig. a) and 3^ (fig. b) 
lattices. 
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Figure 8: The magnetization v as a function of k for H = 0.1. The lattice size is 8^. 
The one-loop and two-loop results for v are represented by the dashed, cf. eq. 
and solid, cf. eq. ( U-i3{ ), curves. 
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Figure 9: The magnetization v as a function of k for several values of k . The five 
curves correspond, from the bottom to the top, to R = 0.1, 0.2, 0.3, 0.4 and 0.5. The 
lattice size is in all cases 8^. The solid lines were obtained by evaluating expression 
for a large number of k, values numerically. 
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is also in this case 8^ and the statistics at each k point the same as in fig. ^ In order 
to better monitor the drop of the magnetization near the FM-FMD phase transition 
we have increased the density of points in that region. The sohd fines represent 



again tfie perturbative result according to eq. ([4.15|) . Tfie agreement between tfie 
numerical data and tfie analytic curve is in all cases excellent. Tfie grapfi sfiows 
tfiat tfie drop to tfie critical point is becoming steeper wfien k is increased. Tfiis 
pfienomenon is a consequence of tfie fact tfiat tfie critical exponent r] in eq. ( |4.10| ) 
decreases witfi increasing k. 

Finally we demonstrate tfiat also tfie volume dependence of tfie magnetization 
data in tfie FM pfiase, cf. fig. is nicely reproduced by tfie perturbative formula 



( f4.15|) . In fig. |iy we fiave plotted f as a function of k for k = 0.2 and five different 



lattices ranging in size from 3'^ to 8^. Tfie statistics at eacfi k value is about 10^ x 
(8/L)2 Metropolis sweeps. Tfie solid lines represent again tfie perturbative results 
according to eq. ( [4.15|) . Tfie two-loop curve agrees nicely witfi tfie numerical data 



down to tfie value of n wfiere tfie analytic curve fias a minimum, but starts to 
deviate wfien k is lowered beyond tfiat value. Tfie two-loop curve increases wfiile 
tfie numerical data continue to fall off. Tfiis sfiows tfiat tfie two-loop formula (|4.15| ) 
is valid only in tfie k interval above tfie minimum. Tfie self-energy S(/c) in eq. (|4.16|) 
diverges at 

1^ = /^min = —2 R {1 — cos27r/L) , (5.8) 

because tfie inverse tree-level propagator Ao^(/c) fias a zero eigenvalue (for some 
nonzero k) at tfiis value of k. Tfiis implies tfiat v in eq. ( [4.15|) approacfies one in tfie 
limit K \ Kjnin- (Lowering k beyond tfiis value would lead to a negative eigenvalue 
of tfie tree-level inverse propagator, and tfiis instability causes q to condense to tfie 
smallest possible value, cf. eq. ( |4.11| ).) Tfie plot sfiows tfiat tfie minimum of tfie two- 



loop curve gets smaller and also narrower wfien tfie size of tfie lattice is increased. 
Tfie minimum drops to zero in tfie infinite volume limit at tfie K-value wfiere a(fi;, n) 
vanisfies, cf. eq. ( |4.19|) . Tfiis k- value coincides witfi tfie one-loop estimate ( |4.20|) 
only in tfie limit k —>■ oo. At k = 0.2 we find tfiat v drops to zero at k ^ 0.03446 
wfiicfi is by about 15% larger tfian tfie one- loop estimate in eq. ([4.20| ). 



From eq. (O) we find tfiat for k = 0.2, Kj^in = -0.6, -0.4, ^ -0.2764, -0.2 



and —0.1 for L = 3, 4, 5, 6 and 8. Fig. |I0| sfiows tfiat tfie two-loop curve approacfies 
one at tfiese values of k,. 



6 Summary and Outlook 

In tfiis paper we fiave calculated tfie pfiase diagram of tfie reduced model for a gauge- 
fixed U(l) lattice gauge tfieory in tfie mean-field approximation. Tfie pfiase diagram 
contains a ferromagnetic (FM), anti-ferromagnetic (FM), paramagnetic (PM) and. 
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Figure 10: The magnetization v versus k for k = 0.2 and several values of L. The 
five curves correspond, from the left to the right, to L = 3, 4, 5, 6 and 8. The five 
solid lines were obtained by computing ( U-IS^ ) numerically on lattices of the same 
size. 
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as a novelty, also a ferromagnetic directional (FMD) phase, where not only the 
U(l) symmetry is broken but also the vector field condenses. The mean-field results 
for the phase diagram were confirmed by numerical simulations. The continuum 
limit of the model corresponds to a continuous phase transition between the FM 
phase and the FMD phase. We have studied the nature of this phase transition 
in lattice perturbation theory and demonstrated that the magnetization vanishes in 
the continuum limit, because of infra-red effects, and that the global U(l) symmetry 
gets restored in this limit. This phenomenon cannot be understood in the mean-field 
approximation. We have shown that the numerical data for the magnetization are 
in good agreement with the results of the weak coupling expansion in the FM phase. 
We have calculated the critical coupling of the FM-FMD phase transition at large 
values of k both in the mean-field approximation and in perturbation theory, and 
find in both cases a small positive value. At k — cx) the mean-field result for the 
FM-FMD phase transition is smaller by about a factor two than the perturbative 
result. 

As a next step we will take up again various proposals for lattice chiral gauge 
theories and investigate whether the problems associated with the strongly fiuctuat- 
ing gauge degrees of freedom can be overcome by gauge fixing. In refs. and [|l^ 
we will show for the case of the reduced abelian Wilson- Yukawa (Smit-Swift) model 
that a) the species doublers decouple in the continuum limit, and b) that the fermion 
spectrum contains only the desired states, namely a massless charged left-handed 
fermion that couples to the gauge field and a massless neutral right-handed fermion 
that decouples. We expect to find similar positive results also for other fermion 
formulations, using a Major ana- Wilson term instead of a Dirac- Wilson term p3| , 
domain wall fermions with waveguide ||^, or staggered fermions |2^]. 



It is challenging to study the U(l) model with gauge fields turned on. It should 
be possible to determine the fermion spectrum in the Coulomb phase and see if it 
remains unaffected at small values of the gauge coupling. A change of the fermion 
spectrum should manifest itself as a new phase transition in the fermion sector. 

It is also important to extend the gauge-fixing approach to nonabelian gauge 
theories. This implies that we first have to specify how to discretize and simulate 
the ghost part of the action ( |2.1| ). The nonabelian case is very interesting because 



in this case we can ask whether confinement emerges at small values of the gauge 
coupling. 
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